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Abstract 

We advance an exact, explicit form for the solutions to the frac¬ 
tional diffusion-advection equation. Numerical analysis of this equa¬ 
tion shows that its solutions resemble power-laws. 
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1 Introduction 


Advection constitutes a transport mechanism of either a substance or a con¬ 
served property by a fluid. This process originates in the fluid’s bulk motion. 
For instance, the transport of silt in a river by bulk water flowing down¬ 
stream. Thermodynamic advected quantities are energy and enthalpy. Any 
substance or conserved, extensive quantity can experience advection by a 
fluid that contains it. Advection does not include transport by molecular 
diffusion. In this work we analyze some features of an evolution equation 
characterizing the combined effects of i) advection and ii) a diffusion process 
described by fractional partial derivatives. The ensuing fractional advection- 
diffusion equation has been recently applied to the study of the transport of 
energetic particles in the interplanetary environment [D [21 El El E], 

There is a considerable body of evidence, from data collected by space¬ 
crafts like Ulysses and Voyager 2, indicating that the transport of energetic 
particles in the turbulent heliospheric medium is super diffusive Bi- Consid¬ 
erable effort has been devoted in recent years to the development of superdif- 
fusive models for the transport of electrons and protons in the heliosphere 
[D EE]. This kind of transport regime exhibits a power-law growth of the 
mean square displacement of the diffusing particles, (Ax 2 ) oc t a , with a > 1 
(see, for instance, 0)- The special case a = 2 is called ballistic transport. 
The limit case a —» 1 corresponds to normal diffusion, described by the 
well-known Gaussian propagator. The energetic particles detected by the 
aforementioned probes are usually associated with violent solar events like 
solar flares. These particles diffuse in the solar wind, which is a turbulent 
environment than can be assumed statistically homogeneous at large enough 
distances from the sun [7]. This implies that the propagator P(x,x',tjt'), 
describing the probability of finding at the space time location ( x , t ) a parti¬ 
cle that has been injected at (x',t'), depends solely on the differences x — x' 
and t — t '. In the superdiffusive regime the propagator P(x, x', t , t ') is not 
Gaussian, and exhibits power-law tails. It arises as solution a non local dif¬ 
fusive process governed by an integral equation that can be recast under 
the guise of a diffusion equation where the well-known Laplacian term is 
replaced by a term involving fractional derivatives [5]. Diffusion equations 
with fractional derivatives have attracted considerable attention recently (see 
PS QIS HH [T21 GL3] and references therein) and have lots of potential applica¬ 
tions ns m- In particular, the observed distributions of solar cosmic ray 
particles are often consistent with power-law tails, suggesting that a superdif- 
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fusive process is at work. 

A proper understanding of the transport of energetic particles in space is 
a vital ingredient for the analysis of various important phenomena, such as 
the propagation of particles from the Sun to our planet or, more generally, 
the acceleration and transport of cosmic rays. The superdiffusion of particles 
in interplanetary turbulent environments is often modelled using asymptotic 
expressions for the pertinent non-Gaussian propagator, which have a limited 
range of validity. A first step towards a more accurate analytical treatment 
of this problem was recently provided by Litvinenko and Effenberger (LE) 
in [6], LE considered solutions of a fractional diffusion-advection equation 
describing the diffusion of particles emitted at a shock front that propagates 
at a constant upstream speed V s h in the solar wind rest frame. The shock 
front is assumed to be planar, leading to an effectively one-dimensional prob¬ 
lem. Each physical quantity depends only on the time t and on the spatial 
coordinate x measured along an axis perpendicular to the shock front. 

In the present contribution we re-visit the fractional diffusion-advection 
equation and provide explicit, exact closed analytical solutions, without ap¬ 
proximations. We also undertake a numerical analysis that shows the these 
solutions resemble power-laws. 


2 Formulation of the Problem 


We deal with the equation 


df d a f 

dt K d\x\ a 


+a % +S{x) ’ 


( 2 . 1 ) 


where t > 0 and f(x,t ) is the distribution function for solar cosmic-rays 
transport. Here the fractional spatial derivative is defined as j[6j 


d a f 

d\x\ a 



T(cx + 1) 


/Qe + O ~ 2 /(» + f(x 

e* + i 



( 2 . 2 ) 


To solve this equation one appeals to the the Green function governed by the 
equation: 


do__ d a g 

dt K d\x\ a 


+ S(x)5(t). 


(2.3) 
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With this Green function, the solution of (12.ip can be expressed as 


t 

f(x,t) = j g(x + at,t)dt. (2.4) 

o 

In this work we obtain the solutions of Eqs. (12.ip and (j2.3[) using distribu¬ 
tions as main tools [16] . 

For our task we use, as a first step, the solution obtained in [6] for the 
Green function through the use of the Fourier Transform given by 

OO 

£(M) = ^ / G{x,t)e~ ikx dx, (2.5) 

—OO 

from which we obtain for G- 

G[k,t) = -n\k\ a G(k,t) + (2.6) 

Z7T 

whose solution is 

£(M) = (2.7) 

Z7T 

where H(t) is the Heaviside’s step function. 


3 Explicit general solution of the equation 

From (12.7p we have for G 

m t) = = m f' ( ; 

2n 2 tt ^ n\ 

n =0 

and, invoking the inverse Fourier transform 

OO 

G(x,t) = ^ [ e~ KWat e ikx dk = 

27T J 

—OO 


Hit) 

2n 


OO 


E 

n=0 


(-1 ) n K n t r ' 


n\ 


z an e ikx dx + 


k an e~ ikx dx 


(3.2) 
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Fortunately, we can find in the classical book of [IB] the results for the two 
integrals of (1.3.21) . We obtain 


&( x , t) = -3- T ( an + 


71—0 


gif(cm+l) g —if(ara+l) 

+ 


_(x + i 0) an+1 (x-i0) an+1 


Using now (12.411 we have for / 


(3.3) 


f(x, t) = / Q{x + at ,t ) dt , 


so that one can write 




(~l) n K r ‘ 


2n *—' n\ 

n =0 


T(an + 1) x 


t 



o 


if (cm+1) 


e 2 


,-i|(an+1) 


+ 


(x + at' + it)) an+l (x + at' — zO) 


an+1 


t n dt 


(3.4) 


According to Eq. (1A.3R of the Appendix, where t > 0, we now obtain for 
/, invoking hypergeometric functions F(an + 1, 2; 3; z) and Beta functions 
B(l,n + 1), 


•«*’*) = ^E 


(-1 ) n K n t 


n nj.n+1 


n =0 


Til 


T(an + 1)£>(1, n + 1) x 


gif (cm+1) 

_(x + z0)" n+1 

g-if (an+1) 

(x - z0)" n+r 


F ( an + 1, n + 1; n + 2; 


at 


F ( an + 1, n + 1; n + 2; 


x + zO 
at 

x — zO 


+ 


(3.5) 


This is the general solution of Eq. 02.ip for the initial condition /(x, 0) = 0. 
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Figure 1: Typical graph for solutions of Eq. (13.51) . for different times t 
(x in arbitrary units). 

4 Results 

Fig. 1 depicts a typical graph for solutions of Eq. (I3.5[) . We do this for 
different times t. The main result of this contribution is that these curves 
are very well approximated by power laws. Figs. 2 (k — 0.01), 3 (k — 0.5), 
and 4 (k = 0.2) show how the above solutions can be adjusted by power laws 
of the form l/x q \ the corresponding lines rotate, for fixed a = 3/2, in the 
plane f(x,t ) vs. l/x q when at changes, as the plots clearly illustrate. In all 
our figures we have taken a = 1 and a = 3/2. 

Figs. 5, 6, and 7 show how the dependence of q with a can, in turn, be ad¬ 
justed in simple fashion. In Figs. 5 and 6, straight lines are employed. In Fig. 
7 we attempt a spline adjustment. It is evident that a weak diffusion-regime 
exists. This regime was conjectured in Ref. [6] by appeal to an approximate 
theoretical treatment of the fractional diffusion-advection equation and here 
amply confirmed by our exact approach. 

What about other a values? The situation does not change. Power-law 
adjustment continues to be possible. Figs. 8 and 9 illustrate this issue. 
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5 Conclusions 


We have provided an explicit analytical solution for and advection-diffusion 
equation involving fractional derivatives in the diffusion term. This equation 
governs the diffusion of particles in the solar wind injected at the front of a 
shock that travels at a constant upstream speed V s h in the solar wind rest 
frame. The shock is assumed to have a planar front, leading to a problem 
with an effective one dimensional geometry, where all the relevant quanti¬ 
ties depend solely on time and on the coordinate x measured along an axis 
perpendicular to the front. 

We obtained the exact solution of the above mentioned equation in the 
^-configuration space (besides the associated formal solution in the A;-space 
related to the previous one via a Fourier transform). 

We conclude from our analysis that the solutions of the fractional diffusion- 
advection equations are essentially power laws, and have numerically found 
a “law” for the behavior of the associated power-exponent q with varying k 
via spline interpolation. 
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Appendix: Some properties of Hypergeometric 
Function 


Using data from D2I we have 


t 


-f-n +1 

dt — -— T — ,. - B(l, n + 1) x 


(x + at' ± ie) an+1 "" (x ± ie) an+1 

at 


F an + 1, re + 1, re + 2; 


x ± ze 


(A.l) 


Then, we have 


lim 


t 


dt = lim 


t n+l 


6^o+ J (x + at' ± ie) an+1 e^o+ {x ± ie) an+l 

at 


B{ 1, n + 1) x 


e—>0+ 

Thus we obtain finally 

r 


lim F an + 1, n + 1, n + 2; 


x ± ie 


(A.2) 


^.n+l 

dt —7— re + 1) x 


(x + at' ± zO) Qn+1 (x ± t0) an+1 

at 


F an + 1, re + 1, re + 2; 


x ± iO 


(A.3) 
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Figure 2: Power-law fitting to the solutions of Eq. (13.5ft for a = 3/2 a = 1, 
and k = 0.01, at different times t 
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Figure 3: Power-law fitting to the solutions of Eq. (13.5ft for a = 3/2 a = 1, 
and k = 0.5, at different times t 
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Figure 4: Power-law fitting to the solutions of Eq. (13. 5 p for a = 3/2, a = 1, 
and k = 0.2, at different times t 
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Figure 5: Possible linear relationship between the power-law exponent q and 

K 
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Figure 6: Possible linear relationship between the power-law exponent q and 

K 
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Figure 7: Spline relationship between the power-law exponent q and k 
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Figure 8: Fractional derivative order a vs. power law exponent q adjusting 
the solutions of the fractional advection-diffusion equation ( k = 0.2) 
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Figure 9: q vs. k for different a—values. 
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